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1. Introduction 


Naturally renewable sources are not constant so their asso- 
ciation with conventional ones permits their uninterrupted 
power generation. Hybrid Energy Systems (HES), combine two 
or more complementary renewable sources like wind and solar 
and one or more conventional sources like diesel generator [1]. 
Generally, most of hybrid systems have a system of energy 
storage [2]. There are many energy storage systems, like, 
electrochemical, inertial and hydrogen [3]. The configuration 
of hybrid system resulting from a design based on the resource 
available and the constraints of uses, in principle the use in 
maximum of the renewable energy resources and the optimiza- 
tion of the power supplied quality. 

Hybrid PV-Wind systems (Fig. 1) offer the most adequate 
solutions for the electrification of remote areas; the combination 
and the ratio of the two types of energy depending greatly on the 
resources locally available in each geographical area. These 
resources can be evaluated only after a period typically one year 
of monitoring of the basic parameters (wind speed, solar radiation 
and temperature), which are necessary for sizing and implementing 
such systems in the respective areas. 


2. Sizing methods and simulation programs review 


The aim for an optimum sizing of stand-alone photovoltaic- 
wind hybrid systems is to guarantee the lowest investment with a 
reasonable and full use of the renewable energies sources, with 
energy storage optimization in terms of load satisfaction. Various 
optimization techniques of sizing have been reported in the 
literature such as proposed by Markvart [4], which studied a PV/ 
Wind hybrid system and determine the sizes of the PV array and 
wind turbine, by using the measured values of solar and wind 
energy at a given location. Another technique for the sizing stand- 
alone hybrid photovoltaic-wind systems has been developed by 
Bagul et al. [5]. This technique has important advantages over the 
original three event techniques: A probabilistic approach is used to 
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get to the results. A general method has been developed to 
determine the sizing and the operation control of hybrid systems. 
The operation control and sizing selection method is based on 
genetic optimization techniques [6]. Kellogg et al. in [7] have 
developed a numerical algorithm for generation unit sizing hybrid 
PV/wind power generating as stand-alone systems. Marwali et al. 
in [8] proposed a methodology for cost calculating production of 
hybrid systems. Also a probabilistic performance assessment 
method was developed by Karaki et al. [9]. A methodology is 
developed to determinate the optimal size of PV-hybrid subsys- 
tems and to optimize the stand-alone system management [10]. 
Diverse methods including probabilistic and/or deterministic 
approaches have been developed to assess the performance of 
these hybrids systems. Energy flow and battery storage system 
was modeled to obtain optimal size of a PV and wind generators. A 
techno-economic analysis for stand-alone PV/wind hybrid energy 
system is presented by Celik [11]. This method is complete by Ai 
et al, which gives more accurate and practical. Also, neural 
network and genetic algorithm may be used and combined for 
sizing and controlling hybrid energy system to giving optimum 
solution [12,13]. 

To evaluate the performance of the hybrid solar/wind systems, 
several software tools are available for designing hybrid systems, 
such as HYBRID2, HYBRIDS, HOMER and HOGA | 14]. 

HYBRID2 is software for simulation hybrid system. It was 
developed by the Renewable Energy Research Laboratory (RERL) 
(University of Massachusetts). HYBRID2 simulations are very 
precise (time intervals from 10 min to 1 h) [15]. 

HOMER (Hybrid Optimization Model for Electric Renewable) 
developed by National Renewable Energy Laboratory. It is a time- 
step simulator using hourly load and environmental data inputs 
for renewable energy system assessment; it facilitates the optimi- 
zation of renewable energy systems based on Net Present Cost for 
a given set of constraints and sensitivity variables. HOMER consists 
of a library including photovoltaic generators, batteries, wind 
turbines, hydraulic turbines, AC generators, fuel cells, electrolyzers, 
hydrogen tanks, AC-DC bidirectional converters, and boilers. The 
loads can be AC, DC, and/or hydrogen loads, as well as thermal 
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Fig. 1. Configuration of the studied stand-alone PV-Wind hybrid system. 
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loads. The user must select the components of the model to 
represent the architecture of his network. For optimization pur- 
poses, technical and financial data for each selected component 
must be entered. It simulates the operation of a system by making 
energy balance calculations for each of the 8760 h in a year. For 
each hour, it compares the electric demand in the hour to the 
energy that the system can supply in that hour, and calculates 
the flows of energy to and from each component of the system. 
The tool also decides for each hour how to operate the generators 
and whether to charge or discharge the batteries (dispatch 
strategies: cycle charging and load following). When sensitivity 
variables are selected as inputs, HOMER repeats the optimization 
process for each sensitivity variable specified. At the end of the 
simulation, the different system configurations are classified by 
their total NPCs [16,17]. 

HYBRIDS, is a Microsoft Excel spreadsheet-based renewable 
energy system assessment application and design tool, requiring 
daily-average load and environmental data estimated for each 
month of the year. It developed at the University of Massachusetts 
and the National Renewable Energy Laboratory. This software uses a 
combined time series/probabilistic method. The time series 
approach is useful for long term predictions but, for an accurate 
diesel dispatch, Hybrid 2 accounts for short term fluctuations in the 
wind and load power. This time steps can vary from 10 min to 
several hours, and here appears the probabilistic or statistical 
approach, applied in each interval. Both the approaches mentioned, 
assure that energy is conserved throughout the entire simulation, 
and that the Hybrid 2 model is internally consistent. The probabil- 
istic approach can be explained due to the fact that the method is 
based on its use of probability density functions for addressing the 
statistics of generation and load short term fluctuations. Using the 
time series/probabilistic method implies that the algorithm 
assumes mean values of the various inputs and internally calculated 
values. Hybrid2 assumes that this values are randomly distributed 
about the mean value used and, furthermore, that the distribution 
values can be described by a probability density function. The 
probability density functions (p.d.f) used are the normal Gaussian 
and in some cases, the Weibull distribution. These p.d.f are used in 
each time step to find expected maximum and minimum values, 
the fraction of time in which those values may be within a certain 
range, and the amount of energy that may be required or available 
in this range [14]. Hybrid2 assumes that the wind speed, wind 
power, and load are all normally distributed over the time step. This 
assumption is based on previous work that showed that the wind 
speed and wind power are approximately normally distributed 
around the mean value over time intervals of approximately 10 min 
[14]. In addition, other work indicates that the electrical load for an 
autonomous diesel grid is also approximately normally distributed 
over a short time interval | 14,15]. Different to HOMER, HYBRIDS can 
only simulate one configuration at a time, and is not designed to 
provide an optimized configuration | 16-18]. 

HOGA (Hybrid Optimization by Genetic Algorithms) is a hybrid 
system optimization program developed by the Electric Engineer- 
ing Department of the University of Zaragoza-Spain. The control 
and the optimization are carried out by means of Genetic Algo- 
rithms, and can be mono-objective or multi-objective. It contains a 
data base of different compounds like photovoltaic panels, wind 
turbines, hydroelectric turbines, fuel cells, as well as batteries, 
regulators, inverters, rectifiers, also a different system loads. HOGA 
Optimization is achieved by minimizing total system costs 
throughout the whole of its useful lifespan, when those costs are 
referred to or updated for the initial investment NPC (Net Present 
Cost). Optimization is therefore financial (mono-objective). How- 
ever, the program allows multi-objective optimization, where 
additional variables may also be minimized: CO, emissions or 
unmet load (energy not served), as selected by the user. Since the 


cost, emissions, or unmet load are mutually counterproductive in 
many cases, more than one solution is offered by the program, 
when multi-objective optimization is required [19,20]. 

The developed software employs a techno-economic approach 
to determine the system that would guarantee a reliable energy 
supply with a lowest investment. This software is mainly char- 
acterized by a simple user interface. Were the user followed five 
clear steps and they can easily enter the specification parameters 
of their system. As it can easily introduce the cost of components 
NPC (Net Present Cost) also costs of maintenance and replace- 
ment (M&R) taking into account the money exchange rate 
variation in the international market. The following table repre- 
sents a resume comparison between different tools under dis- 
cussions (Table 1). 


3. PV-Wind sizing system: Applied methodology 
3.1. Fundamental equations 


For sizing hybrid energy systems some basic equations was 
established [21-23]. 
The first equation describes the photovoltaic power generation. 


Ppy =Npve X Apvg x Gir aí) 


with Apvg(m°) is the surface of PV generator, Npvg the efficiency of 
conversion, given by 


Npvg =r X [1—2 x (Te — T eref)] (2) 


- m Photovoltaic module reference efficiency. 

p Temperature coefficient witch supposed constant for silicon 
solar cells. 

-— Teef Reference solar cell temperature (°C). 

Gir Global solar irradiation (W/m7). 

Te Solar cell temperature it gives by 


NOCT — 20 
Te = Tat (S00) x Gir (3) 


These different parameters are gives by the manufacture of each 
kind of photovoltaic module. 

The second equation describe the wind power generation. 
To determine this power it is necessary to know the wind speed 
[24-26]. So to calculate the wind speed in a specific level a vertical 
extrapolation is establish 
vZ) Z 


E) (4) 


V(Za) Za 


where:v(Z),V(Zq) and a: represent the wind speed in anemometer 
level, the wind speed in the desired level and the coefficient 
characterized wind cutting. 

The Weibull distribution of wind is given by [27] 


r= e- (0 6 


Each wind generator is characterized by a power diagram, which 
depending on functioning speeds. The following equation gives a 
model to calculate the output power of this generator [28,29]. 


Py x gem, Vaem < V < Vn 


Vn — Vdem 
Pw(v) = $ Pn, Vn < V < Vmax (6) 
0, else 
with 


- P,=nominal power, 
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Table 1 
Comparison of commercial programs for sizing hybrid systems. 


Specifications HOMER 
Based method A program based time-step simulation with probabilistic approach and an economic analysis 
Resources Solar, wind, hydro, biomass, fuel 
Load and environmental - Primary, deferrable and thermal load 
data - Hourly load and environmental data inputs 
Compounds - PV modules, wind turbine, hydro, generator, grid, battery, converter, Electrolyzers 
Economic Annual real interest rate, project lifetime, cost of unmet load, system fixed capital cost, system fixed operation and maintenance (O&M) cost, carbon 
tax 
Simulation and analysis - Estimates the cost and determines the feasibility of a system design over the 8760 h in a year 


- Simulates each system configuration and displays list of systems sorted by net present cost (NPC) 
— Performs an optimization for each sensitivity variable 


Optimization Contains the values of each optimization variable that are used to build the set of all possible system configurations 
Results - Simulate various configuration 
- Minimizing total cost referred to the Net Present Cost 


HYBRIDS 
Based method Microsoft Excel spreadsheet-based using combined time series and probabilistic method 
Resources Solar, wind, fuel 
Load and environmental -— Daily-average load 
data - Environmental data estimated for each month of the year 
Compounds PV panels, wind turbine, diesel generators, batteries 
Economic Calculated the NPC 
Simulation and analysis Used 1-h intervals in the simulation 
Optimization No optimization provided 
Results Can only simulate one configuration at a time 
HYBRID2 
Based method A combined probabilistic/time series computer model 
Resources Site parameters as well as time series data of wind, solar irradiation and ambient temperature 
Load and environmental -— Primary, deferrable, optional and heating load 
data - Daily average load and estimated environmental data 
Compounds Wind turbine, PV module, diesel, dump load, battery, converter, synchronous condenser and dispatch strategy 
Economic Costs of the various components as well as economic parameters that are used to evaluate the economic performance of the system 


Simulation and analysis Include values of a number of system variables for each time step. system variables include the power going to each type of load, the unmet load, the 
power produced by each generating unit, the power going into storage, the power conversion losses, the hybrid system diesel fuel consumption, the 
base case system diesel fuel consumption, and the time step energy balance 


Optimization Give the cumulative energy flows and fuel consumption during the simulation run 
Results Simulate only one configuration but the simulations are very precise (time intervals from 10 min to 1 h) 
HOGA 
Based method A program based genetic-algorithm to sizing and control variables of the system add to the probabilities analysis 
Resources Diesel generation, photovoltaic, wind, batteries, Mini-hydro, Fuel cell, electrolyzers, 
Load and environmental - Daily average load 
data - Daily average environmental data 
Compounds PV panels, wind turbines, hydroelectric turbines, fuel cells, H2 tanks, and electrolyzers, batteries, battery charge regulators, inverters, rectifiers and 


AC generators 
Different system loads are possible: electric AC, electric DC, hydrogen, and water pumping load, Thermal load 


Economic Minimizing total cost referred to the Net Present Cost in mono-adjective approach and minimizing total cost referred to NPC and costs of CO and/or 
unmet load in multi-objective approach 
Simulation and analysis — Uses very accurate models for resources, for compounds and for economical calculations 
- Includes an advanced aging batteries model, taking into account corrosion, weighted cycles (effects of SOC, current, acid stratification ...) and 
temperature 


- Includes an advanced model of optimization control strategies 


Optimization — Optimizes the slope of the PV panels, 
- Allows the selection of the currency to be used in the economical calculations 
- Calculates life cycle emissions 
- Allows the buying and selling energy to the electrical grid 


Results More than one configuration is offered, when multi-objective optimization is sought 
— Vaem= Starting speed, The third part of this system is the storage part. The mathe- 
— V,=Nnominal speed, matical model of state of battery charging (SOCpat) is depending 


— Vmax=Maximum speed, on precedent state of charge, and on produced energy by 
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booth wind generator Ew and PV generator Epy, also with load 
demand E; [15,16]. To calculate SOCpat of batteries there is two 
scenarios 


3.1.1. First scenario 
If Ew(t)>E,(t), batteries are in charging process, and the 
instantaneous storage capacity SOC, (tf) is given by: 


SOCpat(6) = SOCpar(t— 1) + (Epy (t) + (Ew(0) — Ex (0) x inv) X Moat D) 
with: 

Epy(t) = Ppy(t) x At 

E,(t) = P,(t) x At 
Ppy (t), Pw (t), Py (t): Are, respectively, PV generator power in an 


instant « t », Wind generator power in the same instant « t >, 
and load demand in this instant « t > [30,31]. 


— At Step time of simulation 

— Niny Inverter efficiency 

— ma Batteries charging efficiency, generally it depends on 
charging current. It is between [0.65 and 0.8]. 


3.1.2. Second scenario 
If Ew(t) < E,(t) we have two possibilities: 


1- If Epy(t) > (EL) —Ew(t))/niny, batteries are in charging state and 
the storage capacity is given by 


E,(t) — Ew(t 
SOCat(t) = SOCyar(t— 1) + (Eo (F0—=)) ta (D 


Ninv 


2- If Epy(t) < (E(t) —Ew(t))/niny, Corresponding to discharging 
batteries, which characterized by the following relationship 
se) ee 1 

Ndisch 


SOCpat(t) = SOCpar(t— 1) + (Eo (10) 


Ninv 


With naison, batteries discharging efficiency, it supposed equal 
to 1. 

In all scenarios the state of batteries charging must satisfy the 
following condition 


SOCpat_ min S SOCpar(O) < SOCpat_ max (11) 


where 


SOCpat_ min ANd SOCpat_ max are the limits states of batteries 
charging. 

SOCpat_ max iS considered like the nominal capacity of the 
storage system, 


The inferior limit is given by 
SOCpat_ min = DOD x Chat_n (12) 


DOD (%) represents the discharging batteries depth. 

The total batteries capacity [W h], depending principle on the 
number of autonomy days ‘Naq, on the produced energy by all 
renewable sources without storage and the daily consumption 
[31,32]. 

We calculate for each day “d” of the year the different between 
the energy demands and produced once. This strategy we can 
permit to identify the variation of SOCpat. 


Eq(t) = E,(t)—Ep(t) (13) 


So 
E(t) = Ew(t) + Epv(t) (14) 
Eq(t) = 24 ,Eq(t) if Ea(Ðmac; sc; 0 (15) 


In the case of energy deficiency, we research the maximal daily 
difference between the energy demands and produced once along 
of a year. To determine the nominal capacity we use the following 
relation [32,33] 


Naa x Max x E,(t) 


Chat n = a 6) 
Ndisch 
The state of batteries charging limits are defines with 
SOCmax = Chat_n (17) 
SOCmin = DOD x Chat_n (18) 


The principle objective of sizing this systems is to determine 
their optimal configuration (size of photovoltaic generator, size 
of wind generator and the capacity of the storage system) to 
satisfy the load demand. Several methods for sizing hybrid 
energy system, was discussed in the introduction. In this 
development the LPSP (Loss of Power Supply Probability) is 
used associated with a technical-economic approach and math- 
ematical models. This method can offered a simulation of 
all configurations possible of the system can satisfy the load 
profile demand. The economic approach is applied to give the 
optimal one. 


3.2. Methodology description 
During the system functioning different cases were distinguished 


a- The load power demand is lower than the produced power by 
the wind generator (P; < Py). In this case the energy surplus 
and the energy produced by the photovoltaic generator will be 
storage on the accumulator park, and the nouvelle state of 
battery charging is calculated like the equation number (07). 
If this state of charge will be higher than the SOCpat max, the 
surplus will be conducting to the auxiliary loads. 

b- The load demand is superior to the wind generator produced 
energy(Pimac; sc; Pw). In this situation the energy deficiency 
(P,— Py) is compensating by photovoltaic generator produced 
energy, and the state of charge is calculated by the Eq. (8). If the 
photovoltaic produced energy cannot satisfy the load, i.e. 
Ppy(t) < (P,(t) — Pw(t))/niny the batteries will discharge to com- 
pensate the energy deficiency with the respect of the dischar- 
ging limit (SOC,,,mac; sc; SOChat. min) and the state of batteries 
charging can be calculated by using the Eq. (9). 


In this case (b), if the batteries cannot satisfy the load demand 
in a desired period « t >», this deficiency is named LPS (loss of 
power supply). It can be calculated with the following relation 


LPS(t) = (PL(t)— Pw) 
x At—(Ppy(t) x At+ Chat(t — 1) — SOCpat_ min) X Ninv (19) 


3.3. Loss of power supply probability (LPSP) description 


This sizing method consists in determining the optimal 
number of the batteries and the photovoltaic modules accord- 
ing to the optimization principle knowing: the reliability, 
which is based on the concept of the probability of loss of 
energy (Loss of Power Supply Probability ‘LPSP’) [35,36], and 
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Data inputs (solar irradiation, temperature, wind speed, Load Profile, : 
LPSP desired, Poy int, Pov sup Pw ints Pw sup, dPpvs dPw, Naa and At Compounds Technical Data 


Day of liberty =1 


P wn E wnt dpw 


Poy n a Pont dppy Pov n= Ppv inf 


Batteries charging 
process, calculus of 
SOCpat (€q7) 


P,y() a (P, (t)-Fy O) Nn 


Batteries charging, Calculus of 
SOCrat (eq8) 


Batteries discharging, Calculus 
of SOC pa (eq9) 


(SOC, 


bat 


Yes No 
> SOC, mx) | 


(SOC,,, < SOC, 


bat — bat_ 


| (SOCpa = SOC. ar max) | | Storage Energie | 


L__, 


(SOG, 4, = SOG a min) 


Calcul de LPS 


No Yes 


Economic 
analysis 


Save all 
configurations 


Optimal 
configuration 


LPSP SLPSPa; 


Fig. 2. Flowchart of the developed algorithm. 


on the cost of the system. ‘LPSP’ is defined as being the fraction charging statement. Otherwise it is defining such as the 
of the deficiency energy and that required by the load. It fraction of all energy losses and the energy demand in a define 
explains the rate of load dissatisfaction, in term of batteries period of functioning (for our study we use t=one year) like it 
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represent on theses expression [22,23] 


F 
LPSP = P,{ Epat(t) < Ebat iii: fort < t} = x t= LPS) (20) 


T PL(® x At 
i.e. the probability that the load statement, at any moment “t”, is 
lower or equal to the minimal edge of the supplied energy in battery 
Ebat-min- Epat(t) is the energy stored in the battery at any moment t, 
expressed in [W h] [23,24]. This method presents the advantage that 
the introducing of the wind generator permits to minimize the cost 
of the photovoltaic stand-alone system, by the minimizing the size of 
the photovoltaic generator and the storage capacities (number of 
batteries) [33,36,37]. 


3.4. The technical-economic developed algorithm analysis 


Based on the LPSP description method, an algorithm was 
developed; 


© Pov infe Ppv_sup and Pw inf, Pw_sup are parameters which repre- 
sent the limits of photovoltaic generator power production and 
wind generator power production, respectively. 

© dPpy and dPw represents the variation steps of photovoltaic and 
wind generator powers. 

e dt represent the simulation step. 


The inputs data of this algorithm are: 


- The average values of solar irradiation, temperature and wind 
speed of each typical day of month in the year. 

- The desired value of the LPSP on one year. 

— Parameters of the different system compounds. 


This algorithm permits the identification of the (Ppy, Pw, Npat) 
triple matrix of all possibilities which can satisfy the LPSP. The 
relationship between these variables is not linear. Only the economic 
analysis can permit to select the optimal sizing of the studied system. 


3.5. The economic analysis and the identification of the optimal 
configuration 


Principally there are three main cost categories: 


3.5.1. The initial cost 
It is related to the purchase cost added to installation cost 


Ci = Pw x Ci w +Ppv x Cipv+Nbat x Ci_pat+Sinv x Ci inv (21) 
with: 


- Ci w: The initial cost of the wind system 

- Ci pv: The initial cost of the photovoltaic system 
— Cj pat: The initial cost of the storage system 

- Cj inv: The initial cost of the inverter 

- Siny: The apparent power of the inverters. 


3.5.2. The maintenance cost 
In the developed software a percentage of the initial cost for 
each compound is gives for a period of a year 


Cm = (Pw x Ci w x Mw+Ppy x Ci pv x Mpyt+Nopat x Cibat 
x Mpat +Sinv X Ci_inv X Minv) x I sys (22) 


- mw: Percentage of the annual maintenance of the wind system 
- Mpy: Percentage of the annual maintenance of the photovoltaic 
system 


— Mpat: Percentage of the annual maintenance of the storage 
system 

— Miny: Percentage of the annual maintenance of the inverters 

- ifgys: Life time of the system [years]. 


3.5.3. Compounds replacing cost 
Each compound of the system have lifetime, it must after this 
lifetime be replaced by another. 


x (If sys i fw) 


Chn= Pw x Ci w Fw +Ppy x Ci pv 
If ys —l Ieys—l 
wef “TF If pv) 4 Noat x one 7 df “F If bat) FSi x Ci inv 
PV bat 
Uf sys ma If inv) 
< Fiw a 


sys» fw lfpv fvat finv: The lifetimes of all the system, the wind 
generator, the photovoltaic generator, storage system and the 
inverter, respectively. Generally the wind generator and the 
photovoltaic generator have lifetime's approaches to the global 
system, so their replaced costs are neglected and can be equal to 
zero [33,34,38]. The global cost of all functioning period is given by 


CSC Cyt, (24) 


In this case, it is very important to introduce the inflation 
ratio, which can take in consideration to the dynamic of money 
changing values. The following figure represents the flowchart of 
the proposed algorithm (Fig. 2). 


4. Sizing developed program: Implementation 


The following section presents the implementation of the 
designed program under windows. These windows described the 
discussed methodology above. This software appears with a 
principal window (Fig. 3) resume all steps for sizing off-grid 


Fig. 3. The five step to sizing PV-Wind hybrid system. 
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PV-Wind systems. These steps are: 
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At each selection of these five buttons, a new window appears; 
allowing the user to define the various parameters characterized 


— Sites settings: introduce the measured environmental data of his system. In the nest section these windows are presented one 


the site of implementation, 


- Load characteristics: establish the load profile, 
- Technical parameters: describe the compounds constitute the 4.1. Sites settings 


by one in order. 


system, 
- Economic parameters: define the parameters of the economic In this window (Fig. 4), the user can select a location from a 
analysis, saved database to implement his system. This database provides 


— System sizing: start simulating and optimization of the system all necessary values for simulation, such air temperature, horizon- 


and give results. 
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Fig. 4. Sites caracterizations. 


Fig. 5. Load profile introducing. 
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4.2. Load characteristics 


This step is crucial; the user specify their power consumption 
profile. Three choices are available: 


- Daily average load: where it should populate the value of his 
office during the 24h of a typical day for each month of the 
year for 12 months. 

- Monthly average load: where it should populate the value of his 
office during a day of each month of the year. 

- Imported data load from an EXCEL file: the user can import data 
profile of a real consumption saved under an EXCEL file. 


The following figure illustrates this (Fig. 5): 

For the daily average load, whenever the user entered a value, 
they must increment time by clicking the “increment” button 
(Fig. 6). 

Under the same window (Fig. 5) the user introduce: the 
minimum and the maximum days of autonomy and the desired 
LPSP. The validation is done by clicking on the “valid” button, 
where the icon red set in green for confirming validation. 


4.3. Technical parameters 


A click on the technical parameters button, a new window 
appears. This window specifically brings all the technical para- 
meters characterizing the system to be sizing. It is unscrewed in 
three fields: 


- Adjustment range production system: the user should specify the 
maximum and the minimum power and the increment step of 
power for each generator (PV and wind), Fig. 7. 

- Parameters of photovoltaic generator: to calculate the power 
generation of PV generator, the user ought to introduce the 
performance of the used panels (commercial modules) especially: 
module efficiency, the nominal operating cell temperature (NOCT), 


Load Consumption Hours 
w 


Fig. 6. Average monthly consumption. 


Variation ranges of production systems 


Minimum Power Maximum Power Increment Step 


SS =a w 
Wind Generator 9000 10100 400 
PV-Generator 2000 3000 100 
Fig. 7. Variation interval of production sources. 
Photovoltaic Generator Parameters 

Nominal Operating Cell Temperature °c 45 
Reference Cell Temperature °C 25 
Temperature Coef Temperature 1/°C -0,0033 
PV module effeciency % 42 


Fig. 8. PV Generator parameters. 


reference or ambient temperature, the temperature coefficient 
B [33], Fig. 8. 

- Parameters of wind generator: under this figure the user should 
specify the ranges of variation of used wind turbines. There are 
three fields to complete, and in each one they inter the power 
range of wind turbine with the speeds characteristics: starting 
speed, Nominal and Maximal speeds (that according to data- 
sheet of the manufacturers) (Fig. 9). 


4.4. Storage settings 


To define the total capacity storage of the system, the user must 
enter the storage capacity in [Ah] of a single battery. Also, all 
parameters described the performance of charging and dischar- 
ging, the depth of discharge and voltage rating must be correctly 
introducing to guaranty a good result of simulation, Fig. 10. 

At the end of a check do not forget to also fill in is the box that 
defines the performance of the inverter use, typically it is around 
90-98% [30], and once the user finish entering all parameters 
specified must be validated by a click on the red icon and the 
green confirms validation. 


4.5. Economic parameters introduction 


In the step number four (see Fig. 3) to sizing these systems, is 
the introducing of the economic parameters which are very 
important to determine optimal configuration of the system 
according to the cost of investment and to the cost of power 
generation. When the user selects this window, they can introduce 
the initial costs of all the compounds constitute the system in $/W 


Wind Generator Parameters 
Starting Nominal Maximum 
Power Curves Parameters Spasd Speed Speed 
w w mjs mjs m/s 
Between 1600 and 4000 45 13 16 
Between 4000 and 7000 45 13 16 
Between 7000 and 9600 45 13 16 


Fig. 9. Wind generator parameters. 


Storage System Parametes 


Storage Capacity [Ah] 200 
Nominal Voltage [v] 12 
Discharging Depth [%] 60 


Batteries Charging Effeciency [%] 80 


Inverter Parameters 


Inverter Effeciency [%] 85 


Fig. 10. Storage system and inverter parameters. 


Valid 
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also the annual maintenance costs and lifetimes of each equip- 
ment of the system, Fig. 11. 


4.6. System sizing: Simulation and optimization 


Arrival in this final and crucial step of “simulation and optimi- 
zation”, the user finishes introducing all parameters characterizing 
their system, a reminder before the design is that the user is 
strongly advised to check the validation of the values recorded 
after each step, verifying that the red icon before the “OK” button 
is set to green in each time. Selecting this button correctly can lead 
the user to the final step (step number five—Fig. 3) to sizing 
correctly the described hybrid system. The following figure repre- 
sents the window that used to starting simulating of the system. 


This program can simulate the operation of a system by making 
energy balance calculations for each hour in a year. For each hour, 
the program compares the produced power in the hour to the 
electric load demand in that hour. In the same time it decides how 
to charging and discharging batteries and whether to operate the 
PV and wind generators. 

The optimization approaches based on LPSP and the economic 
analysis. It becomes after simulating of all configurations possible 
of the system. In the end the program displays of a list of feasible 
configurations, sorted by the optimal one referred to the minimum 
cost. The list can be scanned for other feasible systems. The 
optimum one is put it at the top of the list. 

When the user clicks “start sizing” in the window (Fig. 12), the 
program start the simulation and the optimization, and displays 
the following results described the optimal system configuration: 


Fig. 11. Economic parameters. 


Nominal Power Of Wind Generator (W) 


Nominal Batteries Capacity (Ah) 


Nominal Power Of PY Generator CW) 
[umber of Autonomy Days 


Table ofall = babilities wich's satisfy the LPSP 


9600 8000 


0 5882352941 361,015347255 


0 5882352941 )90,082013922 


9600 10000 | 
9600 12000, 


0 5882352941 319,148680588 


9600, 14000 


0 5882352941 548,215347255 


9600 16000 


o (5882352941 77,282013922 | 


Wind Generator Power Graph 


Number of Batteries Graph 
EEE 


Fig. 12. Screen of simulation results. 
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Table 2 
Developed program specifications. 

Specifications Developed program 
Based method 


Resources 
Load and environmental 


Solar, wind 


A program based on time-step simulation and mathematical models with probabilistic approach and economic analysis 


- Daily average load or Monthly average load or imported data load from an EXCEL file 


data - The average values of environmental data for each typical day of month in the year 


- PV modules, wind turbine, batteries 
- Designs off-grid power systems 


Compounds 


- Simulates each system configuration and displays list of systems sorted by net present cost (NPC), compounds lifetime maintenance and 


Includes an advanced techno-economic strategy to build a set of all system configurations and to minimize total cost referred to the Net Present 


Cost, maintenance and replacement costs with introducing the money exchange variation and inflation rate 


Economic 

replacement (M&R) costs 
Results Simulate various configuration 
Optimization 
Results - Simulate various configuration 


- The optimal power of the wind generator to be used; 
- The optimal power of photovoltaic generator; 

— The optimal batteries capacity; 

- The number of batteries to storage energy; 

- The number of autonomy days; 

- Total net present cost for the whole system. 

— Desired LPSP. 


Other configurations that satisfy the condition specified in the 
value of the LPSP are classified in a table, to fill the calculated 
results from the user presses the button "Add Results", in this way 
the program displays the total list of possible configurations 
preceded by the optimal one. Graphs are provided by clicking 
different buttons: 


- "Draw the graph of PV power generation" 

- "Draw the graph of wind power generation" 
- "Draw the graph of the number of batteries" 
- "Draw the graph of the overall cost". 


The next table gives a resume of the developed program 
specifications, which can be compared with the commercial ones 
(Table 2). 


5. Conclusion 


In this work the LPSP algorithm with a new proposed techno- 
economic algorithm for sizing standalone PV-Wind system was 
implemented. This strategy depend principle on costs study taking 
into account the different compounds costs using in the system, 
their lifetimes also the load profile and the meteorological char- 
acteristics of each installation site. This method was implemented 
under a computer program, which become practical, interactive 
and easy to use. This elaborated program permits to determine the 
optimum size of the battery bank and the PV array for a given load 
and a desired loss of power supply probability taking into account 
the minimum energy cost which is depends generally on invest- 
ment cost, operation and maintenance costs and the depreciation 
period. The competitiveness of the developed program in their 
first version is the simple user interface under windows. 
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